Atomic micromotion and geometric forces in a triaxial magnetic trap 
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Non-adiabatic motion of Bose-Einstein condensates of ru- 
bidium atoms arising from the dynamical nature of a time- 
orbiting-potential (TOP) trap was observed experimentally. 
The orbital micromotion of the condensate in velocity space 
at the frequency of the rotating bias field of the TOP was de- 
tected by a time-of-flight method. A dependence of the equi- 
librium position of the atoms on the sense of rotation of the 
bias field was observed. We have compared our experimental 
findings with numerical simulations. The nonadiabatic fol- 
lowing of the atomic spin in the trap rotating magnetic field 
produces geometric forces acting on the trapped atoms. 

PACS number(s): 03.65.Bz,32.80.Pj,45.50.-j 



Magnetic trapping has become a standard technique 
for confining and evaporatively cooling ultra-cold sam- 
ples of neutral atoms. Up to now, BEC in alkali atoms 
has been achieved in two different kinds of magnetic 
traps: static traps, and dynamic ones, also called time- 
orbiting-potential (TOP) traps [El. In the latter, intro- 
duced in J2|, a quadrupole magnetic field is modified 
by a rotating homogeneous bias field which causes the 
quadrupole field to orbit around the origin, thus elimi- 
nating the zero of the magnetic field at the trap center. 
In the original scheme with cylindrical symmetry p| , the 
bias field rotates in the symmetry plane of the quadrupole 
field, whereas in a later version H the plane of rotation 
contains the quadrupole symmetry axis, which yields a 
completely asymmetric (triaxial) magnetic trap. In the 
present work, precise measurements on the motion of a 
87 Rb condensate in the triaxial TOP trap have allowed us 
to discover new dynamical features, beyond the approx- 
imations previously applied in the theoretical treatment 
of time-dependent magnetic traps. 

In the adiabatic approximation the spins of the mag- 
netically trapped atoms follow the instantaneous direc- 
tion of the magnetic field. Generally speaking, an adia- 
batic approximation is applied when the evolution of one 
set of variables of a system is fast compared to the time 
scale of a second one. While the full dynamics can be 
quite complicated and hard to analyze, in a useful ap- 
proximation the fast motion is solved for fixed values of 
the slow coordinates, and then the slow dynamics is con- 
sidered to be governed by the average values of the fast 
motion. In a classical description, the energy of the fast 
motion calculated for fixed values of the variables of the 
slow system acts as a potential in which the evolution 
of the slow system takes place. In an improved approxi- 



mation, reaction forces act on the slow dynamics. They 
are related to the geometric phase and to the geomet- 
ric gauge fields M. In the classical case, the physical 
origin of the geometric forces has been studied, e.g., for 
atomic systems with spins Hj. For atomic motion in an 
inhomogeneous magnetic field, the slow atomic external 
degree of freedom acts as a "guiding center" for the evo- 
lution of the fast precessing spin variables coupled to the 
magnetic field. The resulting geometric forces are usu- 
ally very small and hard to detect. Analogous terms in 
the calculation of molecular spectra beyond the Born- 
Oppenheimer approximation do, however, lead to mea- 
surable shifts of energy levels, underlining the usefulness 
of the concept Q . 

In the TOP trap, the guiding-center atomic motion 
within the time-dependent inhomogeneous magnetic field 
acquires an additional dynamic feature associated with 
the periodic evolution of the rotating field, feature not 
present in the analyses of ||. Because the frequency of 
the bias field rotation is typically much larger than the 
vibrational frequencies of the atoms in the trap, in an 
additional approximation a time-average over the period 
of the rotating field is usually performed. Application of 
both adiabatic and time-average approximations predicts 
a description of the atomic motion in terms of a time- 
independent potential. A number of theoretical studies 
have considered the role of the nonadiabatic dynamics |j] 
or of the time-dependent potential on the atomic evolu- 
tion in a TOP trap ||7j . Experimental investigations, how- 
ever, have restricted their attention mainly to the atomic 
response to the time-independent adiabatic potential. 

The breakdown of the time-average approximation and 
of the adiabatic spin following approximation, and the 
evidence for additional geometric forces are visible in the 
following measurements on a 87 Rb condensate || in a 
triaxial trap: i) the atomic micromotion at the frequency 
of the bias field, ii) the change of the atomic equilibrium 
position in the trap when inverting the sense of rotation 
of the bias field (|] , iii) the dependence of the rotation- 
sensitive equilibrium position on the applied magnetic 
fields. The dependence of the equilibrium position on the 
sense of rotation of the bias field does not violate any fun- 
damental symmetry law, as the handedness of the Lar- 
mor precession for an atom in a trapped state defines the 
sense of rotation in the atomic evolution. In our experi- 
mental investigation, as for other TOP traps, the typical 
conditions for applying the adiabatic and time-averaged 
approximations are not violated pd[ . A precise detec- 
tion of the micromotion, however, allowed us to perform 
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measurements of the non-adiabatic effects present in the 
TOP trap. We have compared our results with analytical 
models and numerical simulations based on a completely 
classical description, obtaining good agreement. 

Our experimental apparatus |12| is based on a double- 
MOT system with two cells connected by a glass tube to 
capture, cool and transfer atoms into a magnetic TOP 
trap. The latter consists of a pair of quadrupole coils 
oriented along a horizontal axis and two pairs of bias- 
field coils, one incorporated into the quadrupole coils and 
the other along a horizontal axis perpendicular to that of 
the quadrupole coils, as shown in Fig. [I]. After transfer- 
ring the atoms into the TOP-trap, we evaporatively cool 
them first by reducing the radius of the quadrupole orbit 
and finally by applying an RF-field. Once Bose-Einstein 
condensation has occurred in the \F = 2, mp = 2) state, 
we adiabatically changed the quadrupole gradient and 
the bias field to the desired values. The instantaneous 
magnetic field seen by the atoms is: 

B = [2b' x + B Q cosuj t]i + [-b'y + B sinw i]j - b'zk (1) 

where b' is the quadrupole gradient along the z-axis, and 
Bq and wo are the magnitude and the frequency of the 
bias field, respectively. 
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FIG. 1. Schematic of the static and rotating magnetic fields 
confining the condensate, and of the shadow-imaging method. 
Imaging along the j/-axis allowed us to measure the velocity 
of the atoms along the x and z axes. For our trap parameters 
the spatial amplitude of the micromotion in the trap was less 
than 1 /mi, well below our resolution limit of 5 /mi. 

Application of the adiabatic and time-average approx- 
imations leads to an oscillatory motion of the atoms in a 
time-independent potential with a harmonic frequencies 
of 2 : 1 : v% where u> x = b' \f2jjJmBo, 
m being the mass and fi the magnetic moment of an 
atom in the trapped state ||. Lifting the time-average 
approximation yields the atomic micromotion, because 



the atoms experience a time-dependent potential at fre- 
quency loq and its harmonics. Expansion of the potential 
energy U for atoms close to the trap center gives addi- 
tional linear terms AU modulated at the frequency loq 
on top of the effective harmonic potential, with 



AU = pJa (2a;cosa->o^ — ysinujQt) . 



(2) 
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FIG. 2. Dependence of the micromotion velocity vq x on b' 



at wo — 2tt x 10 



in (a) and on wq in (b). In (b), b' 



=329 Gem 1 and Bo =28 G. In (a), for gradients larger than 
200 Gem" 1 , B = 28 G was used, whereas for smaller gradi- 
ents Bo was reduced in order to maintain a small gravitational 
sag. Continuous lines from Eqs. (|^). The inset in (a) shows 
the velocity dependence on the phase of the bias field at which 
the trap was switched off (zero phase corresponds to the bias 
field pointing along the direction of observation) . 

For trap oscillation frequencies well below loq, separa- 
tion of the time scales allows us to write down the atomic 
micromotion at ojq as a small shaking motion on top of 
the secular motion in the trap. Balancing the restoring 
force from the AU potential with the centrifugal force of 
the resulting orbital motion yields a periodic micromo- 
tion with velocity amplitudes 



2/ib' , ub' 
and vq v 



771070 



mu>o 



(3) 



In order to observe the micromotion of the atoms at the 
TOP-frequency, the trap fields were switched off in a time 
(less than 20/Lts) short compared to the period of the 
TOP-field 1 1 3 1 . The instant at which the trapping fields 
were switched off was determined to better than 1% of 
a TOP-cycle. After switching off the quadrupole field at 
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a well-defined phase of the rotating bias field and after 
a variable delay of up to 8 ms, we imaged the absorption 
shadow of the falling cloud onto a CCD-chip. Because 
of the small spatial amplitude of the micromotion, the 
measured position of the atoms along the x-axis after 
the time-of-flight was determined solely by their initial 
velocity vq x in that direction (see Fig. |l|). 

A typical plot of vq x as a function of the bias- field phase 
is shown in the inset of Fig. [| A velocity amplitude 
of 5mms _1 corresponds to a spatial amplitude of the 
micromotion around 100 ran. Also shown in Fig. [2] is 
the dependence of vq x on the quadrupole gradient and 
on the frequency of the rotating bias field. We found 
no dependence of v 0x on B . All of these dependencies 
are in agreement with Eq. (0) (continuous lines in Fig. ||). 
We checked by numerical integration of the equations 
of motion in the adiabatic approximation that the finite 
switching time of the trapping field leads to a reduction 
of the measured micromotion velocity of at most a few 
percent compared to Eq. pL 
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FIG. 3. Centre-of-mass oscillations of atoms in the 
TOP-trap for clockwise and counterclockwise rotation of the 
bias field. The gradient b' and the bias field Bo were 20.8 
Gem -1 and 4 G, respectively. The equilibrium positions for 
the two cases are indicated by the horizontal lines. 

When the sense of rotation of the bias field was re- 
versed, the phase of the measured velocity was shifted 
by 180 degrees, in accordance with the physical picture 
of the atoms being 'dragged' along by the instantaneous 
position of the quadrupole center. On top of this phase 
shift, we observed a shift in the vertical equilibrium po- 
sition of the atoms in the TOP-trap. We achieved a high 
resolution of the equilibrium position by using very small 
condensate clouds and by exciting dipole oscillations of 
the cloud along the z-axis. Typical results of such mea- 
surements for the two senses of rotation of the bias field 
are shown in Fig. [|. Fitting a sinusoidal curve to the 
position data we determined the absolute centre of the 
oscillations to within 5 /^m and the relative position of 
the centers to within 3 /^m. Identical conditions for both 
senses of rotation were ensured by powering the bias coils 



from two phase-locked signal generators and shifting their 
relative phase by 180 degrees to reverse the sense of ro- 
tation whilst keeping all other parameters constant. 
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FIG. 4. (a) shows the dependence on the quadrupole gradi- 
ent b' of the splitting Az oq of the equilibrium positions (filled 
symbols) for opposite senses of rotation of the bias field, and 
dependence of the oscillation frequencies (open symbols). The 
solid and dotted lines are the results of a numerical simula- 
tion without the adiabatic approximation. The dashed line is 
the prediction of Eq.[|; Bo = 4 G for all points. In (b) we plot 
the dependence of the splitting (symbols: data, solid line: nu- 
merical simulation) on the inverse of the bias field strength 
for fixed gradient (b' = 20 Gem -1 ). 

The measured dependence of the spatial separation 
Az cq of the two equilibrium positions on the inverse of 
the quadrupole gradient, 1/6', is shown in Fig. || (a). For 
large gradients, we find a linear dependence of Az cq on 
1/6', whereas for the smallest gradients the splitting data 
deviate significantly from a straight line as 6' approaches 
the minimum gradient, b' in = 15.3 Gem -1 , necessary to 
balance gravity. In Fig. (b) the measured dependence 
of Az eq on the inverse of the bias-field Bq is shown. Also 
plotted in Fig. || (a) is the dependence of ui z on 1/6', al- 
lowing a calibration of 6' with an accuracy of two percent 
over a very large range. 

We have solved numerically the coupled differential 
equations of Jl0[ for the center of mass motion and atomic 
spin within our TOP trap field. The numerical simula- 
tions provided an excellent fit for the results of Fig. ^ 
p4| . To find the equilibrium position, a transient damp- 
ing was introduced in the numerical integration of the 
evolutions for the center-of-mass and/or the spin. We 
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have verified that damping rates corresponding to the 
trap lifetime do not modify the atomic micromotion and 
the non-adiabatic effects. However the role of additional 
damping mechanisms, for instance the viscous heating 
of 0, should be tested. 

To interpret the dependence of the splitting Az oq on 
the field gradient, a transformation to a frame rotat- 
ing synchronously with the bias field may be applied, 
as in |) 10 1 . For the special case of a cylindrical TOP 
trap this transformation completely removes the time de- 
pendence. The frame rotation acting on the spin vari- 
ables introduces a fictitious (geometric) magnetic field 
B g = — ■m-phujo'k/ [j, along the rotation axis, with mph 
the atomic internal angular momentum. In a refined adi- 
abatic approximation the spin is aligned with the effec- 
tive field, i.e. the sum of the real and geometric magnetic 
fields. The occurrence of transverse components of the 
the spin with respect to the real magnetic field lies at the 
heart of the concept of geometric forces. The geometric 
field shifts vertically the plane of the quadrupole orbit. 
The difference Az eq in the equilibrium positions for cw 
and ccw rotations of the bias field is 



Az cq = 2 



Jib' 



(4) 



Gravity shifts the equilibrium position z cq downwards, 
but for a cylindrical trap the above result for Az eq is not 
modified. The expressions for the geometric force of || 
give rise to an equivalent result for the geometric field 
and the splitting Az oq . While for the cylindrical TOP 
trap the transformation to the rotating frame takes care 
completely of the spin dynamics, for the triaxial TOP 
trap the rotating frame transformation does not result in 
a time-independent problem. In fact, at low b' values, 
approaching b' min , the splitting data of Fig. [|(a) deviate 
from the behavior described by Eq. (0). Higher order 
geometric terms derived in |p| lead to corrections several 
orders of magnitude too small to explain these deviations. 
Numerically we found that at zero gravity also for the 
triaxial TOP trap Az cq is well described by Eq. 
This leads us to conjecture that gravity modifies the fast 
correlation between spin orientation and center-of-mass 
motion responsible for the geometric forces. 

In summary, in a TOP trap the micromotion and the 
splitting of the equilibrium positions of the atoms for op- 
positely rotating bias fields are accounted for by a lifting 
of the time-average and adiabatic approximations. The 
dependence of the splitting on the magnetic field gra- 
dient arises from the nonadiabatic characteristics of the 
TOP trap. The micromotion in a TOP-trap is of rele- 
vance for experiments where an accurate control of the 
condensate residual motion on the recoil velocity scale 
(5mms _1 ) is required, e.g. in atom lasers or condensate 
transfer. The observed dependence of the equilibrium 
position on the bias field rotation hopefully motivates 
theoretical efforts to calculate higher order corrections to 



the adiabatic spin evolution in the TOP-trap configura- 
tion. Finally it should be tested, theoretically and exper- 
imentally, whether the condensate response to the par- 
allel transport of the atomic spin in a TOP trap, which 
gives rise to the geometric forces, may influence also the 
internal dynamics of the condensate. 
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The switch-off time of the bias field, limited by the qual- 
ity factor of tank circuits was around 150 fis. However, 
fast switching of the bias field was not required as the 
atoms are no longer trapped once the quadrupole field is 
extinguished. 

We should mention that the agreement between simula- 
tions and experiments becomes less satisfactory for very 
elliptical bias field in the x — y plane. 
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